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ABSTRACT 

We investigate the properties of the two-dimensional distribution and kinematics of 
ionised gas and stars in the central kiloparsecs of a matched sample of nearby active 
(Seyfert) and inactive galaxies, using the SAURDN Integral Field Unit on the William 
Herschel Telescope. The ionised gas distributions show a range of low excitation re- 
gions such as star formation rings in Seyferts and inactive galaxies, and high excitation 
regions related to photoionisation by the AGN. The stellar kinematics of all galaxies in 
the sample show regular rotation patterns typical of disc-like systems, with kinematic 
axes which are well aligned with those derived from the outer photometry and which 
provide a reliable representation of the galactic line of nodes. After removal of the non- 
gravitational components due to e.g. AGN-driven outflows, the ionised gas kinematics 
in both the Seyfert and inactive galaxies are also dominated by rotation with global 
alignment between stars and gas in most galaxies. This result is consistent with previ- 
ous findings from photometric studies that the large-scale light distribution of Seyfert 
hosts are similar to inactive hosts. However, fully exploiting the two-dimensional na- 
ture of our spectroscopic data, deviations from axisymmetric rotation in the gaseous 
velocity fields are identified that suggest the gaseous kinematics are more disturbed 
at small radii in the Seyfert galaxies compared with the inactive galaxies, providing a 
tentative link between nuclear gaseous streaming and nuclear activity. 

Key words: galaxies: active - galaxies: spiral - galaxies: Seyfert - galaxies: kinematics 
and dynamics - galaxies: structure. 
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1 INTRODUCTION 

It is now widely accepted that central supermassive black 
holes (SMBH) are ubiquitous in bulge-dominate d galaxies 
IIFerrarese fe Merrittl l2000l : iGebhardt et~ai] l2000h . but on- 
going nuclear activity is observed only in 10 to 20 per 
cent of nearby gala xies, 3 per cent at the Seyfert level 
l|Kewlev et al.l 120061 ) suggesting black hole re-ignition via 
accretion of material in the central parts of the galaxy is 
required. A key unanswered question is whether the igni- 
tion mechanism is related to the galaxy host properties in 
particular the fuel transportation mechanism. Spiral galax- 
ies contain a rich supply of potential fuel and the search 
for a mechanism to remove angular momentum and drive 
material towards the nucleus has been the focus of many 
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studies. Numerical simulations predict efficient gas radial 
inflows due to non-axisymmetric perturbations in the grav- 
itational potentialsuch as bars or driven by galaxy inter- 
actions (|Wadall2004l h which may help to bring gas close to 
the galactic nucleus. However, most of past imaging stud- 
ies failed to find any significant differences between galaxies 
hosting an Active Galactic Nucleus (AGN) or not, on spa- 
tial scales that encompass nearby companions/galactic in- 
terac tions (|De Robertis et al.|[l998l ; ISchmittll200ll: iHo et al 



2003), stellar bars or nuclea r spirals (Mulchac y fe Regan 
1997: Laurikai nen et al.l 12004 iMartini et alJl2o6a T In con 



Hunt fe Malkanl feo04l ) conducted a NICMOS/HST 



trast, 

imaging study of the circum-nuclear regions of 250 nearby 
galaxies and found a significant excess of isophotal twists in 
the circumnuclear regions of Seyfert 2 galaxies, compared to 
other categories (HH/Starburst, Seyfert 1, LINER, inactive 
galaxies). This result may hint at an evolutionary scenario 
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to explain the different states in which a nearby galaxy is 
observed, and could suggest the presence of identifiable dy- 
namical differences between Seyfert and inactive galaxies in 
the central kpc regions. 

Probing the dynamics of Seyfert galaxies requires 
spectroscopic data. Single aperture and long-slit spec- 
troscopy studies are clearly inadequate to investigate 
the complex structures observed in the central kpc of 
Seyfert galaxies. Two-dimensional spectroscopy (integral- 
field spectroscopy, IFS hereafter) is therefore a pre- 
requisite to study the dynamics of the gaseous and stel- 
lar components. Early work with such instruments pro- 
vided two-dimensional gas and stellar kinematic maps 
of the central parts of a few near by active galaxies, in 
optic al wavelength : e.g. NGC4151 ([M cdiavilla fe Arribasl 
Il995l). NGC3227 iMediavilla fe Arribasl Il993r i. NGC 1068 
i Garcia-Lorenzo et all Il997l ) with IN TEGRAL (an IFS 
based on optical fibers), or NGC 2110 [jPsrruit et al.ll2004l ) 
using OASIS (an IFS using a microlens array), as well 
as in the near infrared (NIR): The Circinus galaxy ob- 
served wi th the NACO spectrograph and SINFONI a t 
the VLT dPrieto et al.ll2004l; iMueller Sanchez et al.l l2006h . 
NGC 3227 l|Davies et al.l I2006T ) observed with SINFONI. 
Multiple gaseous systems and kinematic perturbations are 
sometimes revealed, but the corresponding field-of-view 
(FOV) and/or the angular resolution were often too small 
to disentangle the AGN-related and galactic disc line emis- 
sions. Statistical conclusions cannot be drawn from these 
detailed studies of individual Seyfert galaxies, originally tar- 
geted for their complex nuclear properties and lacking any 
control inactive galaxy comparison. More recenly, the molec- 
ular gas in the central part of a larger (though still small) 
sample of Seyfert galaxies was observed using the Plateau 
de Bure Interferometer in the course of the NUGA sur- 
vey (|Garcia-Burillo et al.ll20o"3T ). revealing the potential role 
of gravity torques in the feeding of the inner 100 pc (e.g. 
iGarcfa-Burillo et alj|2005h . Gas reponds non-linearly to de- 
viations from axisymmetry so a more direct link with the 
gravitational potential is still difficult to establish, due to the 
lack of large-scale stellar kinematic maps for these galaxies. 
Six nearby Seyferts were also observed with GMOS / Gemini, 
enabling the mapping of the ionis ed gas and stellar kin emat- 
ics in the central few arcseconds (jBarbosa et alj|2006r i. IFS 
such as SINFONI /VLT rea ched regions even closer to the 
nucleus l|Davies et al.|[2006l ). thanks to the use of adaptive 
optics. Again, such observations are useful for studying the 
kinematics in the nuclear regions, identifying the putative 
molecular torus around the SMBH, retrieving the character- 
istics of the central engine such as the mass of the SMBH, or 
studying the environment of the active nucleus. However, ob- 
served kinematics at these relatively small scales are rather 
difficult to interpret in the context of the host galaxy, consid- 
ering the significant influence of non-gravitationally driven 
processes (e.g. outflows). 

The SAURON IFS, mounted on the William Herschel Tele- 
scope (La Palma, Spain) has FOV (33" x 41" in low reso- 
lution mode) that is large enough to observe ionised gas 
under the influence of the galaxy host potential, and a high 
enough spatial sampling to still probe the AGN-related emis- 
sion in the inner parts. These data are then well suited to 
study both the stellar and gaseous kinematics in the above- 
mentioned context of fueling. A number of key studies were 



recently performed with SAURON pro viding unique datasets 
for a sample of nea r by early-type |de Zeeuw et all |2002| ; 
lEmsellem et al.l [20041: ISarzi et al.l l2006l) and sp iral galaxies 



( Ganda et al. 20061 : iFalcon-Barroso et al.l 120061 ) . As a spin 
off, the distribution and kinematics of the gas and stellar 
components were also obt ained and studied fo r some well 
known Seyf e rt (N GC 1068. lEmsellem et al.ll2006l , NGC 5448, 
iFathi et all 120051) or spiral inactive galaxies like M 100 
jAllard et al.ll200o ). These works provided evidence for the 
existence of ionised gaseous inward streaming in the inner 
few kpc, hinting at a mechanism for transporting gas in the 
circumnuclear regions and further. However, no survey has 
been yet pursued to map the two-dimensional kinematic and 
morphology in the central parts of a well-selected sample 
of Seyfert and inactive galaxies, in order to search for and 
quantify potential dynamical differences between these two 
populations at these spatial scales. 

We therefore designed a new IFS survey for a well- 
defined sample of active and control inactive galaxies, se- 
lected from the RSA sample with nuclear classification of 
iHo et all l| 19971 ). This program aims at probing the galactic 
potential of the host galaxies by comparing the properties 
of both the ionised gas and stellar components in the cir- 
cumnuclear regions of Seyfert galaxies. We observed these 
galaxies with the SAURON IFS, which provides a spatial cov- 
erage sufficient to map regions from the central 33" x 41" ~ 
few kpc down to the inner 1" ~ 50 pc at 10 Mpc. In 
Sect. we present the sample selection, the observations 
and data reduction. In Sect. [3l we present the associated 
stellar and ionised gas maps. A kinematic analysis is pre- 
sented in Sect. U We then briefly discuss our results in a 
more general context and conclude in Sect. [5] 



2 OBSERVATIONS AND ANALYSIS 
2.1 Sample Selection 

The observations presented here form part of a larger mul- 
tiwavelength campaign to investigate host galaxy struc- 
ture and kinematics on size scales ever-closer to the nu- 
cleus, aimed at identifyi ng or eliminating possible triggering 
and f ueling mechanisms l|Mundelj|2o'ol : lGarci'a-Burillo et all 
2005). Our master sample consists of 39 Seyfert galaxies se- 
lected from the RSA catalogue paired with 39 control galax- 
ies with carefully matched optical properties: Bt, Hubble 
type, inclination and, where possible, V ays (^ 4000 kms -1 ). 
3D spectroscopic imaging with the VLA and WSRT is un- 
derway to map the large-scale neutral hydrogen Hi distribu - 
tion and kinematics for these galaxies (|Mundell et al.ll2007l ). 
Current Hi imaging interferometers however cannot rou- 
tinely resolve structures smaller than ~ 20" (i.e. r < 700 pc 
for our sample); studying the distribution and kinematics 
of ionised gas and stars with the integral-field spectrograph 
SAURON provides the missing link between the large-scale Hi 
disc properties and the very central nuclear regions. 

For our SAURON study, we selected a distance-limited 
sub-sample of Seyfert+control galaxies with V sys < 
1600 kms" 1 to ensure that Fe stellar absorption lines lie 
well within the spectral band for the full SAURON field-of- 
view. Measurement of these Fe lines is critical in the pres- 
ence of emission lines from ionised gas that contaminate Mg 
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Figure 1. R-band Digital Sky Survey images of the sample galaxies. The '(S)' or '(C)' on the right of the object names stand for Seyfert 
or Control galaxy, respectively. Each Seyfert galaxy is displayed on the left of its associated control galaxy, except for NGC 1068 and 
NGC3227 (first two panels) for which no control has been observed. The orientation is such that North is up and East is left, the bar 
located at the bottom right corner of each panel corresponds to the spatial length of 5 kpc. Overplotted on each image is the position of 
the SAURDN field of view for that galaxy. 



lines (see Sect. [2]3}- In total, our sub-sample comprises 15 
pairs of Seyferts+inactive galaxies. We completed observa- 
tions of 7 pairs and two well-known Seyferts (NGC 1068 and 
NGC 3227) with no control galaxy data. NG C 1068 has been 
studied in detail by Emscllc m et al.l (|2006h and is included 
here for completeness, while weather constraints prevented 
observations of the control galaxy for NGC 3227. Table [1] 
lists the properties of the sample galaxies. The V-band ab- 
solute magnitudes of the Seyfert nuclei in this sub-sample 
span the full magnitude range of the 25 brightest Seyfert 
nuclei in the RSA catalogue (-20.9 < My < -16.0), thus 
offering a representative selection of Seyfert activity. 



2.2 SAUR0N Observation and Data Reduction 

Observations of our sample were carried out between 2001 
and 2004 using the integral-field spectrograph SAUR0N at 
the 4.2m William Herschel Telescope (WHT) at La Palma, 
Spain. Table [2] summarizes the exposure time for each point- 
ing and lists the equivalent seeing of final reduced data. 

The low spatial resolution mode of SAUR0N was used, 
providing a field-of-view (FOV) of 33" x 41" with a square 
sampling of 0.94 arcsec per spatial element (lens). This de- 
livers about 1500 spectra simultaneously per object which 
cover the spectral range 4825-5380 A, with a resolution of 
4.2 A. This wavelength range includes a number of impor- 
tant stellar absorption lines (H/3, Fe5015, Mg&, Fe5270) and 
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Table 1. Properties of our sa mple. (1) Pairs identifier, (2) Galaxy name, (3) Hub ble Type (NED), (4) Numerical morphological type 
(LEDA), (5) Distance i n Mpc llHo et al.lll997}) , (6) D 25 in arcmin llHo et alj |l9970. (7) disc inclination in degrees, (8) Total apparent 
magnitude Bt in mae |Ho et al.lll997h . (9) Systemic vel ocity in kms -1 (NED), (10) V-band absolute magnitude of the Seyfert nuclei 
My in mag llHo et alj|l997fl , fill Spectral classification dHo et al.lll997Tl: H=HII nu cleus, S=Seyfert nucleus, L=LINER, T=transition 
object, 'T2' implies t hat n o b road Ha emission line s w as detecte d dHo et al.lll997h. { 12) References for the disc inclination values. 1: 
iGarcia-Gomez et al.l J2002h 2: iMundell et alj lll995h, 3: [Ho et alj Jl997h . 4: iKassin et al.l ( |2006h . 5: Ijogee et al] (l2002f) , 6: iPedlar et ahl 
jl992h . 7: lErwinl l|2005l l. 8: lTullvl dl974h . 9: lRozas et alJ 



ionised gas emission lines (H/3, [O ill] and [N i]). More details 
on th e SAUR0N spectrograph can be found in iBacon et al.1 
|200ll ). Fig. [T] presents the SAURDN field-of-view overlaid on 
R-band Digital Sky Survey (DSS) images of our galaxies. 

The data of the 15 galaxies were reduced using the 
dedicated XSAUR0N software and an automatic pipeline 
available within the SAURON consortium l)Bacon et al.l l200ll : 
Ide Zeeuw eraLll2002f ). The main steps include: bias and dark 
subtraction, extraction of the spectra using a mask, wave- 
length calibration, low-frequency flatfielding, cosmic rays re- 
moval, homogenisation of the spectral resolution in the field, 
sky subtraction (using 146 sky spectra 1.9 arcmin away from 
the main field). Th en the flux calibrat i on wa s applied as ex- 
plained in detail in iKuntschner ct al. (2006). The flux cali- 
brated individual exposures are then accurately centred with 
respect to each others and merged. The final merged dat- 
acube is sampled onto a squared grid with 0'.'8 x 0'.'8 pixels. 

Finally, the point spread function (PSF) of each merged 
exposure was determined by comparing the reconstructed 
SAURON intensity distribution with HST/WFPC2 images. 
The SAURON PSF is modeled by a single 2D Gaussian, whose 
parameters is determined by minimizing the differences be- 
tween the HST/WFPC2 images, convolved by this Gaussian 
and the SAURON reconstructed images. The derived seeing 
values (Full Width at Half Maximum) are listed in Table [2] 

2.3 Derivation of Stellar Kinematics 

In order to ensure the reliability of the stellar kinematic mea- 
surements, the merged datacubes were spatially binned us- 



ing th e Voronoi 2D binning algorithm of Ca ppellari fc Copinl 
(2003), creating bins with a minimum signal-to-noise ra- 
tio (S/N) = 60 per bin. The stellar contribution to the 
SAURON spectra is potentially contaminated by emission 
lines: H/3A4861, [O m] AA4959, 5007 and [N i]AA5198,5200, 
therefore we first identified and masked out the spectral re- 
gions that are significantly contaminated by emission. The 
stellar kinematics are measured on each spectrum in our 
binned datacubes usin g the penalized pixel fit t ing (p PXF) 
method developed by ICappellari fc Emselleml ([2004) . The 
line-of-sight velocity distribution (LOSVD hereafter) is pa- 
rameterised by Gauss-Hermite series and the algorithm de- 
termines the mean velocity V, the velocity dispersion a and 
the higher Gauss-Hermite moments /13 and hi which mini- 
mize the differences between the observed spectrum and a 
stellar template spectrum convolved with the corresponding 
LOSVD. A low-order polynomial is included to account for 
small differences between the galaxy and the template spec- 
tra. A first estimate of the kinematic parameters is done 
using a single star spectrum as template. Then using this 
initial estimate, an optimal stellar template is derived for 
each individual s pectrum, via t he use of a large library of 
stellar templates l|Vazdekislll9"99h . We finally iterate by mea- 
suring best-fitting values at each position for the velocity V, 
the velocity dispersion a, hz and hi using this time the opti- 
mal templates obtained from the previous step. The stellar 
kinematic parameters were determined using this method 
automatically implemented in the SAURON data reduction 
pipeline for our whole sample except for the Seyfert 1 galax- 
ies: NGC 3227, NGC 4051 and NGC 4151. For these galaxies, 
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Pairs Name dates T exp FWHM PSF 

NGC (s) (arcsec) 

(1) (2) (3) (4) (5) 

1068 2002 January 300 + 3x1800 1.9±0.1 

3227 2004 March 3x1800 1.4±0.2 



1 2655 2004 March 6x1800 1.4±0.1 

4459 2001 March 4x1800 1.6±0.1 



2 3627 2004 March 4x1800 1.6±0.1 

5806 2004 March 2x1800 2.1±0.1 



3 4051 2004 March 6x1800 1.4±0.2 

5248 2004 March 5x1800 1.9±0.1 



4 4151 2004 March 3x1800 2.0±0.3 

2985 2004 March 3x1800 1.2±0.1 



5 4579 2004 March 4x1800 1.9±0.1 

3351 2004 March 4x1800 0.9±0.1 



6 5194 2004 March 4x1800 1.2±0.1 

5055 2004 March 6x1800 2.4±0.1 



7 6951 2003 August 3x1800 1.4±0.1 

5248 2004 March 5x1800 1.9±0.1 



Table 2. Details of the exposures of our sample. (1) Pairs identi- 
fier, (2) NGC number, (3) date of observation (4) Exposure time 
in sec, (5) Seeing, Full Width at Half Maximum in arcseconds. 

the automatic pipeline alone is inadequate due to the pres- 
ence of a broad H/3 emission line in the inner few arcseconds 
(see Fig. [2]c). The spectral regions that are masked have to 
be carefully defined, and then we fit the stellar kinematics 
of these 3 galaxies interactively. 

The AGN continuum is not separated from the stel- 
lar component using this technique. The stellar kinematic 
derivation of the Seyfert 2s is a priori not affected by the 
AGN continuum since the central engine is hidden from di- 
rect view. The AGN continuum would however affect the 
derivation of the velocity dispersion of the Seyfert Is in the 
central regions where it is dominant. These regions have 
been excluded from the stellar kinematic analysis and sub- 
sequent discussion. 

2.4 Derivation of Gaseous Distribution and 
Kinematics 

The stellar continuum resulting from the stellar kinematic 
fit was subtracted from the original data, providing pure 
emission line datacubes. The wavelength range of our ob- 
servations includes the H/3A4861, [O ill] AA4959, 5007 and 
[N l]AA5198,5200 emission lines. The parameters of these 
emission lines (intensity, mean velocity and FWHM) were 
derived from Gaussian fittin g using the fit/spec software 
developed bv lRoussetl (|l992T ). 

The fit was performed on two systems of emission lines : 
the first consisting of the H/3 and [N i] lines, and the sec- 
ond, of the [O in] lines. The H/3 and [O in] lines were fitted 
independently in order to detect differences in the kinemat- 
ics of the two lines. Since the [Nl] lines are significantly 
weaker than the H/3 and [O ill] lines, their kinematics can 
not be constrained independently, therefore this doublet is 
fitted together with the H/3 line system. Within each sys- 



tem, the lines were assumed to share the same velocity and 
FWHM. Constraints were applied on the parameters: the 
line ratio [O m]A4959/[0 m]A5007 was assumed throughout 
to be equal to 2.88, and 0.5 < [N i] A5198/ [N i] A5200 < 1.5. 
Moreover the velocity V and the velocity dispersion FWHM 
of the lines were bounded: V around the systemic velocity 
of the galaxy, and the FWHM by the spectral resolution 
of SAUR0N FWHM SAUR0N = 4.2 A. Usually the line profiles 
were simple enough to be fitted automatically. A first fit was 
done using one single Gaussian profile and this automated 
fit has been visually controlled for each galaxy. An example 
of the emission lines fitting results is given in Fig. [5]a. A 
small number of galaxies show complex emission lines pro- 
files and required a manual fit with additional emission line 
components. Additional components for [O ill] lines, H/3 line 
or both were required in the case of five active galaxies: 
NGC 3227 (Fig. gib), NGC 4051, NGC 4151 and NGC 5194. 
These additional components are certainly associated with 
the Narrow Line Regions (NLR) of these galaxies. An unre- 
solved central H/3 broad component has also been added in 
the case of the Seyfert 1 galaxies (NGC 3227, NGC 4051 and 
NGC 4151) which corresponds to the Broad Line Regions 
(BLRs). The FWHM of this component is 2600 kms -1 , 
1500 kms" 1 and 3100 kms" 1 for NGC 3227, NGC 4051 and 
NGC 4151, respectively, and the spatial extent of the ra- 
dial profile of the BLRs (Fig. [3} is consistent with the PSF 
derived above, within the derived uncertainties (Table [2}. 
Finally, in the case of NGC 4051, Fell blends are detected 
as shown in Fig.[2]d, and these emission lines were therefore 
added in the fit for this galaxy. 



3 RESULTS 

In this Section, we present the distribution and kinematic 
maps of the stellar and ionised gas components for our sam- 
ple of galaxies. Figures [3a] to I4hl present our SAUR0N maps: 
stellar continuum, [O ill] and H/3 intensity distributions and 
emission line ratio [O m]/H/3, as well as stellar and gaseous 
kinematics (velocity and velocity dispersion). In the case of 
the Seyfert 1 galaxies, the line ratio [O ill] /H/3 is computed 
using the narrow H/3 component. The galaxies are displayed 
by pair: the Seyfert on the top, its control galaxy below, ex- 
cept for NGC 1068 and NGC 3227 which are shown together 
(Fig. I4a[) . All the maps are oriented so that the outer pho- 
tometric major-axis of the galaxy is on the horizontal axis. 
To display the gas maps, we show only the reliable emission, 
i.e. when the ratio of the fitted amplitude to the surrounding 
noise is larger than 3. The [Nl] doublet is detected in most 
of our galaxies, but it is very weak, so it is not discussed fur- 
ther. The gaseous kinematic maps (velocity and dispersion) 
correspond to the [O ill] emission lines kinematics, except for 
the inactive galaxy NGC 4459 (Fig. [4b]) for which the H/3 line 
was used since it is slightly more extended and has a better 
signal-to- noise ratio (Fig. [3c]). The H/3 emission line veloc- 
ity and velocity dispersion maps are shown in Appendix [B] 
([Oiii] corresponding maps for NGC 4459). 

Finally, NGC 4459 and NGC 6951 were o bserved as part 
of other programmes (|de Zeeuw et al.l 12002') with multiple 
fields to provide a mosaic (Fig. [TJ. For consistency with the 
other galaxies in our sample, we extracted images corre- 
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wavelength (A) wavelength (A) 



Figure 2. Examples of spectra and their corresponding fits for 4 galaxies, (a): NGC 2985 (inactive), (b): NGC 3227 (Seyfert 1.5), (c): 
NGC4151 (Seyfert 1.5), (d): NGC 4051 (Seyfert 1.2). In each panel, the black line on the top shows the galaxy spectrum and the red 
line, its corresponding fit which is composed of the best-fitting stellar template added to the emission line fit. Below this spectrum, the 
pure emission line spectrum is shown in black. The corresponding emission lines fit is overplotted in red. The spectra were taken in the 
central pixel of the FOV, except for NGC 3227, for which the spectrum was taken at 4" North-East from the center. 
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Figure 3. Radial normalised profiles of the BLR for NGC 3227, NGC 4051 and NGC 4151 (filled blue circles). The red solid line 
corresponds to the best-fit seeing value for each of the three galaxies (see Table |2j, the two red dashed curves correspond to the 
associated upper and lower limits at 2a. 
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sponding to one single SAURON field of view (FOV) from the 
mosaiced exposures for these two galaxies. 

In the following we present the general properties of our 
sample. Detailed descriptions of the maps for each galaxy 
can be found in Appendix [XI 

3.1 Stellar and Ionised Gas Distribution 

For each galaxy, the stellar continuum maps (Figs. I4al4hl 
first panel for each galaxy) were derived by integrating 
over the full wavelength window the spectra correspond- 
ing to the optimal stellar template obtained as explained 
in Sect. 12.31 We constructed the ionised gas intensity maps 
and the [O m]/H/3 lines ratio maps directly from the fit of 
the emission lines spectra. 

3.1.1 Stellar Continuum Distribution 

Half of the galaxies in our sample present symmetric stel- 
lar continuum maps with regular isophotes. NGC 2985, 
NGC4151 (Fig. gejl and NGC 4459 (Fig. E9 show rather 
round central features while the isophotes are more flattened 
for NGC 3227 (Fig. Ejj}, NGC 4051 (Fig. H), NGC 5806 
(Fig. ic) and NGC 6951 (Fig. EE}. The latter four galaxies 
host a large scale stellar bar, their stellar continuum reflect- 
ing the elongation of the bar in the central regions. 

Five Seyfert 2 galaxies and one inactive galaxy present 
more complex stellar distribution maps. Isophotal twists 
or irregular isophotes are observed in five Seyfert 2 galax- 
ies: NGC 1068 (Fig. ga), NGC2655 (Fig. EE}, NGC3627 
(Fig. go), NGC 4579 (Fig. El and NGC 5194 (Fig.gjf. The 
inactive galaxy NGC 5055 presents an asymmetric structure 
(Fig. |4g[ ). The surface brightness is higher on the North- 
West side of the FOV, the emission line flux being absorbed 
by the dust on the South-East (see Fig. [TJ. The predom- 
inance of irregular isophotes and twists in the circumnu- 
clear regions of Seyfert 2 galaxies compared to Seyfert 1 
or inactive galaxies h as been described quantitatively by 
iHunt fc Mailanl (|2004l ). 

Finally, two inactive galaxies, NGC 3351 and NGC 5248 
exhibit a circumnuclear ring at a radius of ~ 5" which 
corresponds to 550 pc for NGC 3351 and 200 pc for 
NGC 5248, respectively. These struc tures have been de- 
scribed bv lDevereux et ail l| 19921 ) and lLaine et ail (|200ll i 

3.1.2 Ionised Gas Distribution 

Ionised gas is detected in all the galaxies of our sample, 
over the full SAURON FOV except for NGC 4459 (Fig. lib)) , 
where emission lines are very weak outside the inner 10". 
In this Section we describe the maps of the gas properties 
mainly focusing on regions which are not dominated by non- 
gravitational motions driven by the active nucleus. 

A variety of structures can be seen in the [O ill] and H/3 
intensity maps: 

• Spiral-like structures are seen in two Seyfert galax- 
ies: NGC 1068 (Fig. EH H/3 intensity map) and NGC 4579 
(Fig. El ■ Such structures corr espond well with t he nu clear 
molecular spirals ob s erved by ISchinnerer et al] (|2000f ) and 
iGarci'a-Burillo et ail (|2005h in NGC 1068 and NGC 4579, re- 
spectively. 



• Circumnuclear rings are found in three inactive galax- 
ies (NGC 3351 Fig. E3 NGC 5248 Fig. gd] and NGC 5806 
Fig. ic} and the Seyfert 2 NGC 6951 (Fig. EE}, which cor- 
respond with known ring-like star for ming structures in the 
central regions of these four galaxies l|Devereux et al.lll992l : 
lLaine et al1l200ll : ICarollo et aill2002l . respectively). 

• In five other galaxies, asymmetric structures are ob- 
served. The emission line distribution of NGC 4151 is elon- 
gated from the centre to the South- West side of the FOV 
(Fig. l4e} in agreemen t with the high-exci tation emission line 
feature described by IPerez et ail (| 19891 ) . corresponding to 
the Extended Narrow Line Regions (ENLR) of this galaxy. 
The Seyfert 2 galaxy NGC 2655 (Fig. EH} presents a hot-spot 
15" away from the centre to the East and a lane 10" West of 
the centre elongated along the Sout h/North direction consis- 
tent with the polar ring observed bv lKeel fc Hummell |l988). 
Asymmetric and irregular structures observed in the ionised 
gas intensity maps of NGC 3627 (Fig. ic} or NGC 5055 (gg} 
are certainly due to the presence of dust in the circumnuclear 
regions. In the case of NGC 5194 both the stellar contin- 
uum and the H/3 emission line intensity p resent irregularities 
consi stent with the dusty nuclear spiral (|Peeples fc Martini 
2006), while [O ill] emission lines traces an outflow structure 
associated with the AGN. 

• Finally three galaxies (two inactive NGC 2985, Fig. Hel 
NGC 4459, Fig. ES and the Seyfert 1 NGC 4051, Fig. Ed} 
show regular round gaseous distributions. 

3.1.3 Ionised Gas Line Ratios 

As expected, the Seyfert galaxies present higher [O m]/H/3 
ratios than their associated inactive galaxies in the central 
kpc regions. The maximum measured ratio for the inactive 
galaxies is of the order of 2, while for Seyfert galaxies it can 
reach values up to 10 or 20. For the Seyfert galaxies, the 
emission line ratio reaches its highest value in the central few 
arcseconds associated with the active nucleus. Some regions 
away from the central engine are also characterised by high 
[O iii]/H/3 values (^> 1) corresponding to high excitation re- 
gions, such as in the ionisation cone of NGC 1068 (Fig. I4a[) . 
the bubble of NGC 5194 (Fig. |4g[ ) or the elongated gaseous 
feature of NGC 4151 (Fig. Ee}. Finally the ring-like struc- 
tures seen in the gaseous maps of NGC 3351, NGC 5248 and 
NGC 5806 show low emission lines ratios ([Om]/H/3 < 0.1) 
consistent with star formation. 

3.2 Stellar and Ionised Gas Kinematics 

Outside the regions dominated by non-gravitational motions 
associated with AGN-driven outflows, all of our galaxies 
show stellar and gas velocity fields with a global regular ro- 
tation pattern (second panel of each row in Figs. l4alto l4h|) . 
More complex dynamical structures are observed in some 
cases, especially in the ionised gas velocity fields. 

3.2.1 Stellar Kinematics 

Most of the galaxies in our sample present a regular stel- 
lar velocity field, the orientation of their kinematic major- 
axis being constant over the FOV and their minor-axis per- 
pendicular to it (e.g. NGC 2985 Fig. El NGC 4459 Fig. EB 
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NGC5055 Fig. |4g[ ). A few active galaxies show some de- 
partures from axisymmetry e.g. the S-shaped zero-velocity 
line in NGC 1068 (Fig. I4a|l and twisted kinematic major- 
axis in NGC 3627 (Fig. [4c}. In all of our sample galaxies, 
the stellar kinematic major-axis is rather well aligned with 
the outer disc photometric major-axis (see Sect. 14. lTTj) . The 
observed stellar kinematics within the SAURON field of view 
are also clearly dominated, in all cases, by rotational mo- 
tion, as evaluated via a global measurement of V/a (Dumas 
et al., in preparation), confirming the disc-like nature of the 
galaxies in our sample. 

For some of our galaxies, the stellar ve locity dispersion 
decre a ses inwards. Th e se so-c alled <r-drops (Emscllc m et al.l 
l200ll ; IMarauez et"ail 120031 ) are observed in at least 
five Seyfert galaxies NGC 1068, NGC 2655, NGC 3227, 
NGC 4051, NGC 4151 and NGC 6951 and one inactive 
galaxy NGC 5248. These central cr-drops have been found 
to be common in SO to Sb spira l galaxies (lEmsellem et alJ 
|2004| ; iFakon-Barroso et al.ll2006l ; IMarauez et al.ll2003l) and 
recently iGanda et al.l (|2006l ) showed that such structures 
are present in galaxies as late as Sd. Central velocity dis- 
persion drops are thought to be associated with dynami- 
cally cold structure, like discs and star formation regions 
|Wozniak et al.ll2003h . 

3.2.2 Ionised Gas Kinematics 

The ionised gas velocity maps are dominated by rotation 
and for all of the galaxies the gas rotates faster than the 
stars (Figs. l4"al to l4hl second panels of the second rows). Since 
stars follow collisionless orbits, their mean velocities are in- 
deed ex£e£ted_to_l2e_l_ower_than the local circular velocity 
(see iBinnev fc Tremaind Il987h . The gas velocity fields are 
more distorted and present richer structures than the stellar 
ones. Evidence for significant deviations from axial symme- 
try exist in almost all the galaxies: radial change of orien- 
tation of the major-axis (NGC 3227, Fig. I4ap . wiggles along 
the kinematic minor-axis (NGC 2655, Fig. gb] , NGC 4579, 
Fig- 12}, S-shaped feature (NGC 1068, Fig.gaJ, or more com- 
plicated structures (e.g. NGC 3627, Fig. [2c] or NGC 5194, 
Fig. [4g] ) . For all of our galaxies, the H/3 and [O ill] velocity 
fields present very similar structure, though [O ill] velocity 
fields have generally better signal-to-noise ratio than the H/3 
ones, except for NGC 4459 (Figs. l4aland lBla[) . We also ob- 
serve highly misaligned kinematic major-axes of the ionised 
gas with respect to the stellar ones in some of our galaxies 
such as NGC 2655 (Fig. HQ and NGC 4579 (Fig. Hi}. These 
misalignments will be analysed in the following Section. 

The gaseous velocity dispersion maps present no 
particular structure for all of the inactive galaxies 
and two Seyferts (NGC 3627, Fig. gc] and NGC 4051, 
Fig. gd|). For eight Seyfert galaxies (NGC 1068, NGC 2655, 
NGC 3227, NGC 3627, NGC 4151, NGC 4579, NGC 5194 
and NGC 6951) and three inactive galaxies (NGC 2985, 
NGC 4459 and NGC 5055) a values increase inwards. In the 
inner 5", the dispersion can reach values 200 km s _1 higher 
than in the outer part of the FOV as in NGC 1068 (Fig. RE)) . 
NGC 2655 (Fig. |3b]) and NGC 4579 (Fig. |4f|) also present 
high dispersion values (~ 230 km s" 1 ) associated with the 
gas distribution structures. For NGC 4151, gaseous disper- 
sion velocity maps exhibit high values in the inner arcsec- 
onds (about 280 km s _1 ) corresponding to the central en- 



gine. Outside the very inner part and inside the ENLR, the 
velocity dispersion shows lower values (between 50 and 130 
km s _1 ). Then, a increases in the outside parts of the map 
for this galaxy, correspo nding to the location of a dusty ring 
jVila-Vilaro et al.lll995l ). 



4 COMPARISON OF STELLAR AND 
GASEOUS KINEMATICS 

In the previous Section we described the SAURON maps of 
the stellar and ionised gas components for the Seyfert and 
inactive galaxies. For all of our galaxies, the stellar velocity 
fields are dominated by rotation, showing regular isovelocity 
contours and their photometric and kinematic major-axes 
seem globally aligned. In comparison with the stellar maps, 
however, the gaseous distributions are more complex and 
the kinematics more perturbed with this effect appearing 
more pronounced for the Seyferts compared to their con- 
trol galaxies. The identification of this qualitative difference 
between gaseous and stellar distributions and kinematics in 
the inner kpc of active and inactive galaxies has important 
implications for the triggering and nuclear fueling. In this 
Section, we therefore provide a more quantitative analysis 
of the stellar and gaseous velocity fields and the photometric 
major- axis of the galaxies. We first compare the major-axis 
orientation of the stellar velocity field to that of the outer 
galactic disc. Then we quantify the differences between stars 
and gas in Seyfert and inactive galaxies as seen in our SAURON 
maps. We compare the global orientation of the stellar and 
ionised gas kinematic maps to probe any significant differ- 
ences between the two components. Then we analyse the 
kinematic maps in more detail, by computing the kinematic 
parameters (systemic velocity, dynamical centre, PA, rota- 
tional velocity) as a function of radius, using a simple rep- 
resentation of a two-dimensional thin disc in rotation. 



4.1 Global Kinematic Major-Axis Orientation 

The stellar and gaseous velocity fields can be first simply 
parameterised by their major-axes. Misalignments between 
the kinematic major-axis of the stellar and ionised gas com- 
ponents may reveal the presence of non-axisymmetric struc- 
tures or decoupled components such as central discs. Here, 
we compare such misalignments in our sample of Seyfert and 
control galaxies, to characterise global kinematic differences 
in their central kpc. 

The kinematic axis orientations were determined for 
both the stellar and ionised gas velocity fields. The global 
position angles (PAs) of the kinematic major-axis were ob- 
tained by minimizing the differences between our SAURON 
velocity fields (stars and [O ill]) and a bi-antisymmetric 
representation of thes e fields. This method is detailed by 
iKrainovic et alJ (|2006T ) in their Appendix C and we used a 
specific implementation written by Michele Cappellari. The 
values for the measured kinematics PAs and misalignments 
are listed in Table [3] 
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Figure 4a. SAURQN maps for the Seyfert galaxies NGC 1068 (top panels) and NGC3227 (bottom panels). For each galaxy the first row 
shows (from left to right) the reconstructed continuum map, the stellar velocity field, the stellar velocity dispersion map and H/3 intensity 
distribution. The second row consists of the [O III] intensity velocity and velocity dispersion maps and the line ratio [O III] /H/3 map. All 
the intensity maps are in units of 10~ 16 erg s _1 cm -2 . The stellar and emission lines intensity maps are overlaid on the stellar and 
ionised gas kinematics. Velocities and velocity dispersions are in kms -1 . The color scale is shown on the right hand side and the cut 
levels are indicated in the top right corner of each map. All the maps are orientated so that the photometric major-axis of the galaxy is 
on the x-axis. The long and short arrows on the right of each galaxy name show the North and the East directions, respectively. 
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Figure 4b. SAURON maps for NGC 2655 (top panels) and its control NGC4459 (bottom panels). See caption of Fig.|4a]for details. 
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Figure 4c. SAURON maps for NGC 3627 (top panels) and its control NGC 5806 (bottom panels). See caption of Fig. [4a] for details. 
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Figure 4d. SAURON maps for NGC 4051 (top panels) and its control NGC 5248 (bottom panels). See caption of Fig.|4a]for details. 
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Figure 4e. SAURON maps for NGC4151 (top panels) and its control NGC 2985 (bottom panels). See caption of Fig.|4a]for details. 
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Figure 4f. SAURON maps for NGC 4579 (top panels) and its control NGC 3351 (bottom panels). See caption of Fig. [4a] for details. 
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Figure 4g. SAURON maps for NGC 5194 (top panels) and its control NGC 5055 (bottom panels). See caption of Fig.|4a]for details. 
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Figure 4h. SAURON maps for NGC 6951 (top panels) and its control NGC 5248 (bottom panels). See caption of Fig. [4a] for details. 
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(blue circles). The dashed black horizontal line corresponds to a 
difference of 10° . 



4-1.1 Stellar Kinematic and the Orientations of the Line 
of Nodes 

In the case of an axisymmetric mass distribution, the direc- 
tion of the stellar kinematic major-axis should coincide with 
the line-of-nodes (LON) of the galaxy, while in the case of 
a triaxial potential the kinematic and photometric major- 
axes can depart from each other due to projection effects. 
We first therefore check if the stellar kinematic major-axis 
in the central regions and the outer photometric major-axis 
are aligned. The photometric major-axis position angle is 
computed from ellipse fitting on the R-band DSS images 
for each of our galaxies. These derived values are then com- 
pared with published ones. For most of our galaxies, the 
PAs values found in the literature are in good agreement 
with those derived from the ellipse fit, within the error bars, 
in which case we took our measured values as the pho- 
tometric PAs. However, a few galaxies present significant 
discrepancies between our fitted values and the published 
ones. This is the case for NGC 1068 and NGC4151 which 
present a weak oval bar and a faint, almost circular outer 
disc, respectively. The measurement of the photometric PA 
of NGC 5194 is uncertain given that this galaxy is nearly 
face-on (i=20°). We therefore used reliable values de rived 
from Hi kinematics for NGC 4151 jPedlar et alJll992h and 
optical kinematics for NGC 5194 i|Tullvlll974T )~as the pho- 
tometric PAs. Table [3] lists the resultant photometric PA 
values for the sample galaxies. The differences between the 
global stellar kinematic and the photometric PAs are plotted 
in Fig. HI 

Our first qualitative impression is confirmed: the stel- 
lar kinematic PAs are parallel to the photometric ones 
within 10° for all of our galaxies except for NGC 5194 
and NGC 1068 which present a difference of about 20° and 



15°, respectively. In the case of NGC 1068, the stellar kine- 
matic major-axis PA corresponds to the average PA of 
the kinematic axis as f itted on the large-scale Hi data by 
iBland-Hawthorn et al.l |l997), while the outer photometri c 
major-axis lies at a PA of about 80° (|Emsellem et alj|2006h . 
Besides the two specific cases of NGC 5194 and NGC 1068, 
for all of our galaxies the major-axis orientation of the cen- 
tral part of the stellar velocity field is therefore a reliable 
measurement for the line-of-nodes of the galaxy. 



4-1.2 Global Kinematic Misalignments between Stars and 
Gas 

Fig. [6] presents the global gaseous versus stellar kinematic 
PAs and the distribution of kinematic misalignments be- 
tween the stellar and gaseous components in our sample. 
Most galaxies in our sample follow the simple one-to-one 
relation (left panel). Five Seyfert galaxies show an abso- 
lute misalignment larger than 20°, and for all the inac- 
tive galaxies, the misalignment between gaseous and stel- 
lar kinematic major-axes is less than 20°. Among the 5 
Seyfert galaxies with strong kinematic misalignments, two 
are Seyfert 1 galaxies (NGC 3227 and NGC 4051) for which 
the BLR is detected and one is NGC 5194 which presents 
an outflow structure in the gaseous maps on the North- 
ern side (Fig |4g[ ). The other two galaxies are NGC 2655 
and NGC 4579. NGC 2655 is known to host off-planar gas 
in the central regions (|Erwin fc Sparkell2003l ) and an inner 
molecular spiral is present in the nuclear region of NGC 4579 
jGarcia-Burillo et alj|200"5l ), a spiral-like structure is also ob- 
served in our ionised gas distribution maps (Fig.|4f|. In order 
to take out any influence of the AGN itself on the result of 
our analysis, we masked the BLR for the Seyfert 1 galax- 
ies and the prominent outflow region of NGC 5194 in the 
stellar and gaseous velocity fields. Then we determined new 
values for the kinematic position angles of the stars and 
gas. These values of the kinematic misalignment are listed 
in Table [3] for these four galaxies, and we show them as 
non filled red diamonds in the right panel of Fig. HJ] These 
new values are different from the PAs determined without 
masking the regions where the dynamics is dominated by 
the AGN. For NGC 4151 and NGC 5194, the kinematic mis- 
alignment between the stellar and ionised components is 
smaller after masking such regions. However, it is the op- 
posite effect for NGC 4051: the difference between the stel- 
lar and gaseous kinematic major-axis orientation is larger 
after masking the BLR. For NGC 3227, there is no signif- 
icant difference between the two values. After this second 
analysis, the same five Seyfert galaxies have still kinematic 
misalignments larger than 20°. 

Although about half of our Seyferts sample presents 
strong misalignments between stellar and gaseous kinemat- 
ics while none of the inactive galaxies shows such differ- 
ences, no statistically-significant difference is found between 
Seyfert and control-galaxies. A Mann- Whitney U test shows 
that the two distributions are identical. By symmetrizing the 
velocity fields to look for global differences in disc orienta- 
tion, this zeroth order analysis shows that our sample size 
is too small to be able to separate statistically the two dis- 
tributions of Seyfert and inactive galaxies, despite the fact 
that half of the Seyferts present kinematic misalignments 
larger than 20°. 
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3 


115±2 


130±1 


5 ±4 


15±3 


4 


4151 


SI 


22 


5 


19±3 


43±1 


1±3 


24±4 




2985 


inactive 


-3±1 


3 


-2±2 


-3±1 


1±3 


1±3 


5 


4579 


S2 


96±2 


3 


99±1 


154±1 


3±5 


55±2 




3351 


inactive 


-168±2 


3 


-164±1 


-154±1 


6±3 


10±2 


6 


5194 


S2 


-190 


6 


-165±5 


-134±1 


25±2 


30±6 




5055 


inactive 


102±2 


3 


99±1 


106±1 


3±4 


7±2 


7 


6951 


S2 


137 


7 


145±2 


142±1 


7±3 


3±3 




5248 


inactive 


110±2 


3 


115±2 


130±1 


5±4 


15±3 



Table 3. Results of the kinematic analysis of the stellar and gaseous components. (1): Pair identifier; (2): NGC number; (3): Activity; 
(4): Photometric position angle (PA) in degre es. The values for NG C 41 51 and NGC 5 1 94 co rrespond to the line of nodes PA. (5): 
O rigin of the values for the p hoto metric PA: 1 : Emscllc m et al. I <2006ft . 2: iMundell et al l l|l995 j), 3: R-Band DSS images ellipse fitting, 
4: IChristopoulou et all jl997l) . 5: iPedlar et alj dl992T l, 6: iTullvl dl974T ), 7: iRozas et alj d2002h : (6): Stellar Kinematic position angle 
(PA), in degrees; (7): Gas Kinematic position angle (PA), in degrees; (8): Difference between the photometric and the stellar kinematic 
PAs, in degrees APA p / lot _j.j„ = \PA p i lot —PA s tars\; (9): Difference between the stars and gas kinematic PAs, in degrees APAkj n = 
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Figure 6. Kinematic misalignments between stars and gas. Left panel: global gaseous versus stellar kinematic PAs. The one-to-one 
relation is shown as a black line, active galaxies and inactive galaxies are shown as filled red diamonds and filled blue circles respectively. 
Right panel: distribution of differences of global kinematics PAs of stars and gas PA s tars — PA gas for each pair of galaxies, Seyfert 
(red diamonds) and inactive galaxies (blue circles). The non-filled red diamonds correspond to the PA derived by masking the BLR, of 
the Seyfert 1 galaxies NGC 3227, NGC 4051 (pair 3) and NGC 4151 (pair 4) and the outflow bubble region for NGC 5194 (pair 6). The 
dashed black horizontal line corresponds to a kinematic misalignment of 20° . 
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4.2 Radial Variations of Kinematic Properties 

In the above analysis, global parameters were extracted from 
the two-dimensional data, and therefore details of the ob- 
served kinematic structures and small-scales variations of 
the PA are not accounted for. Since we have the full two- 
dimensional dynamical information, we can analyse the stel- 
lar and ionised gas velocity fields quantitatively over the 
SAURON FOV and then study the kinematic properties as 
a function of radius. Fourier expansion of the velocity field 
along ellipses is a powerful technique to unveil the kinematic 
perturbations and t herefore to give information on the grav- 
itational potential (|Schoenmakers et al.l Il997l ; IWong et al.l 
|2004 iKrainovic et al.ll2006h . This formalism will be used in 
a second paper (Dumas et al. in preparation) on our SAURON 
data associated with large scale radio data obtained with the 
Very Large Array (VLA). For the current paper, we estimate 
the kinematic parameters of our stellar and gaseous veloc- 
ity fields by modeling the galactic disc as a two-dimensional 
thin disc, in pure circular rotation with no vertical velocities. 
The line-of-sight velocities can then be written as: 

Vios(R, = Vsys + Kot(-R)cos$sini 

+ V;. ad (7?)sin$sini (1) 

where V sys is the systemic velocity of the galaxy, V ro t 
and Vrad are the rotational and radial velocities in the galac- 
tic plane. They depend only on the radius R. The disc incli- 
nation i ranges from 0° (face on) to 90° (edge on). Finally, 
the polar coordinates in the plane of the galaxy (R, $) are 
related to the observable coordinates {x, y) on the sky plane 
by: 

_________ _ _________ 

(2) 

sin$ = -0 - Xc) cos $q - (y - Y c ) sin $ 
Rcosi 

(X C ,Y C ) are the coordinates of the dynamic centre and 
$o is the PA of the projected major-axis of the disc measured 
with respect to North (counter-clockwise). 

We use the expression of the velocity (Eq. [1} with the 
tilted-ring method in order to derive the stellar and gaseous 
kinematic parameters. Our technique is i nspired from th e 
ROTCUR routine in the GIPSY package (lBegemanlll987l ). 
The galactic disc is divided into concentric ellipses deter- 
mined by 6 parameters : the centre (X C ,Y C ), position angle 
$0, inclination i, the offset velocity V sys and the rotational 
velocity V ro t- For this first iteration, V ra d is set to 0. The 
width of each ring is chosen to be large enough to have a 
sufficiently high number of pixels so the fit is good, and 
small enough to assume the rotational velocity is constant 
along the ring. The width of the rings is finally chosen to be 
equal to the seeing. All these parameters (except V ra d) are 
adjusted iteratively for each ring, by fitting the expression 
of the line-of-sight velocity of Eq.[T]to the observed velocity 
field: 

(i) Initial values for V ays , &o and i were taken from previ- 
ous published studies or the NASA Extragalactic Database 
(NED). We fixed the dynamical centre to coincide with the 
photometric centre, or with the location of the central radio 
continuum peak for the Seyfert galaxies. 



(ii) A first fit is done, by fixing $0 and i while allowing 
V S ys, (X c , Y c ) and V ro t to vary as a function of R. The mean 
values of V sys , X c and Y c were taken as improved estimates 
for the systemic velocity and the dynamical center of the 
galaxy. 

(iii) Then, we fit the position angle $0 and V ro t with V sya 
and (X C ,Y C ) fixed at the values derived from the previous 
step, and the disc inclination fixed at its initial guess 

These three steps were followed for the stellar velocity 
fields of our sample. For the gaseous velocity fields, we im- 
ported the values of V 3ya and (X c , Y c ) derived from the first 
step of the fit applied on the stellar components. Then we 
fitted the PA and the rotational velocity as for the stars. For 
the Seyfert 1 galaxies and for NGC 1068, the central few arc- 
seconds were excluded since they are highly contaminated 
with emission from the active nucleus (BLR or NLR). 

The results of this analysis are presented in Figs. I7al 
and 1 7b I The stellar and gaseous kinematic PAs are plotted 
as a function of radius. The values of the global stellar and 
gaseous kinematic PAs derived in Sec. 14.11 are indicated on 
each plot as a dotted red and dashed blue lines, respectively. 
All angles are given relative to the photometric major-axis 
PA. These plots reveal the change of orientation of the stellar 
and gaseous kinematic major-axis orientation as a function 
of radius. The global values of the kinematic PAs derived 
in the Sec 14.11 correspond only to a small part of the FOV 
on which the kinematic orientation is globally constant. For 
NGC 3627, the global stellar and gaseous kinematic PAs de- 
rived in the previous Section correspond to values outside 
R > 15" (Fig. [Taj, and for NGC 5806, our first method is 
not sensitive to the variations of PAs between 10 and 15" 
(Fig. [Taj. In the case of NGC 2655, the kinematic PA of the 
ionised gas rises abruptly for r < 10" from the value found 
by our global method to a value close to the stellar kinematic 
PA. 

In order to quantify the variations in the radial profiles 
of the kinematic PAs, we have used the difference between 
the maximum and minimum values of the kinematic PAs. 
For each galaxy, the quantity APA=PA maa; -PA m i n is com- 
puted for the gaseous and stellar components. To compare 
the Seyfert hosts to their associated inactive galaxies, the 
maximum and minimum values of the kinematic PAs were 
derived on regions corresponding to the same spatial length 
for the two galaxies in each pair. If the velocity fields of the 
Seyfert and its associated inactive galaxies extend outside 
r = 1.5 kpc, APA is derived inside this radius. If one or both 
of the paired galaxies extend less that 1.5 kpc, the minimum 
radius value of the two galaxies is used. The central few arc- 
seconds where residual AGN-induced perturbations may be 
present are excluded. The resulting APA values for both the 
gaseous and stellar components, and for all galaxies in our 
sample (Seyfert and inactive) are provided in Fig. [8] Also 
shown is the average APA of each component for the inac- 
tive galaxies (blue lines) and the Seyfert galaxies (red lines) . 
For the latter, the average is made excluding NGC 1068 and 
NGC 3227 since these two Seyfert do not have associated 
control galaxies. 

The first result of this analysis comes from a compar- 
ison between the gaseous and stellar components: for both 
the Seyfert and inactive sub-samples, APA is almost always 
higher for the gas than for the stars. While the stellar kine- 
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Figure 7a. Kinematic position angles of the stars (red diamonds) and of the gas (blue circles) components as a function of radius. The 
first row shows the Seyfert galaxies NGC 1068, and NGC 3227. Each subsequent row presents a pair of galaxies: the Seyfert on the left, 
and the inactive galaxy on the right. The regions related to the BLR of the Seyfert 1 galaxies (NGC 3227, NGC 4051 and NGC 4151) 
and to the NLR of NGC 1068 have been excluded from the fit. The dotted red and dashed blue horizontal lines represent the values of 
the global position angle of the kinematic major-axis for the stars and gas respectively. These values were derived from the velocity field 
symmctrization method (see Section l4.1l l. The values are relative to the photometric PA (see Tabic [3J. 
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Figure 7b. Fig. I7al continued. 
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Figure 8. Distribution of the amplitude of variations of the stellar (left panel) and the gaseous (right panel) kinematic PAs in the inner 
1.5 kpc (see text for details). Active galaxies and inactive galaxies are shown as filled red diamonds and filled blue circles respectively. 
The dotted blue and dashed red lines correspond to the average of PA max -PA m i„ for the inactive and Seyfert galaxies respectively, for 
the stellar (left panel) or gaseous (right panel) components. NGC 1068 and NGC3227 have been excluded to compute the average values 
corresponding to the Seyfert galaxies, since they do not have an associated inactive galaxy. 



matic PAs radial variations are restrained within about 10° 
and 20° for the inactive and Seyfert galaxies, respectively, 
such variations of the gaseous PAs reach more than 40° in 
the case of four Seyfert galaxies (NGC 1068, NGC 2655 and 
NGC 3627 and NGC 4051). The average APA of the gaseous 
component is about ~ 3 times greater than the average APA 
of the stars for the Seyfert galaxies as well as for the inactive 
galaxies. The second result is obtained when comparing the 
Seyfert and inactive galaxies: for both gas and stars, APA 
is almost always larger for the Seyferts than for their con- 
trol inactive galaxies. The average APA for the Seyferts is 
about 2.5 times greater than the average APA for the in- 
active galaxies, considering either the stellar or the ionised 
gas components. 



5 DISCUSSION AND CONCLUSION 

5.1 Accretion Rates and Kinematic 
Misalignments 

In this section we highlight the relationship between kine- 
matic disturbances in the galaxy and the accretion rate of 
the nuclear black hole. Th roughout this p aper we use the nu- 
clear type classification of lrlo et al.l (|l997T ). Figure[9]presents 
our galaxies in the diagnostic diagram Log([0 m] A5007/IL3) 
vs Log([N n] A6583/Ha), as well as the bo undaries betwee n 
the different excitation classes adopted bv lHo et all l|l997h . 
The separation between Seyfert and LINERS objects is not 
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Figure 9. Diagnostic diagram for the galaxies of our sample, 
Log([Om]A5007/H/3) vs Log([N n]A6583/H»). The emission line 
ratios are from Ho et al. (1997). Active galaxies are shown as filled 
red diamonds and inactive galaxies as filled blue circles. Finally, 
the classification boundaries from Ho et al. (1997) are shown as 
black dashed lines. 
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as sharp as depicted by this diagnostic diagram. We can no- 
tice that NGC 5055 is close to be a Seyfert while NGC 3627 
and NGC 4579 are Seyferts and close to the LINERS and 
Transition objects. The activity classification of the sam- 
ple nuclei therefore does not necessarily accurately repre- 
sent the nuclear black hole accretion rate. We have then 
estimated the Eddington rates and the corresponding mass 
accretion rates of the individual nuclei from four relatively 
independent methods. The results are summarized in Ta- 
ble |U The Black Hole masses in this Table (column 3) are 
from direct dynamic measurements (e.g. NGC 1068), from 
reverberation mapping results (e.g. NGC 4051, NGC 4151), 
or estimated from the central stellar velocity dispersion in 
our SAURON maps via t he well-studied M-ct relationship 
(e.g. iTremaine et~aDl2002r i . We have used four observational 
estimates for the Eddington ratio of the nuclei. The first is 
the optical continuum luminosity of the galaxy nucleus from 
HST images (column 4 of Table[4| , interpret ed as optical syn- 
chrot ron radiation from the base of the jet (|Chiaberge et al.l 
120051 ). No bolometric correction is applied so the result 
should be considered a lower limit to the Eddington ratio. 
Of course, contamination of the nuclear optical continuum 
from star formation could cause the AGN-related luminos- 
ity value to be lower. We have also used the hard (2-10 
keV) X-ray luminos ity, the nuclear [O m] luminosity (from 
iPanessa et al]|2006l and lHo et ailll997T) and the nuclear ra- 
dio power to estimate the bolometric luminosity of the nu- 
cleus (columns 5 to 7 of Table [4]). In the case of the [Oiii] 
luminosity, we tested the estimator used f or active galax- 



ies in the SDSS 



3500 x L 



|q jjjj (|Heckman et al.l 

120041 ') . This formula resulted in too hi gh values of Lboi for our 
sample. The factor 3500 derived bv lHeckman et all (|2004| ) 
was based on the empirically derived ratios Lsooo/Ljq m ] 
= 320, where L5000 is the monochromatic c ontinuum lumi- 
nosity APa at 5000Aand Lboi/Lsooo — 10.9. Hcc kman et al.l 
(2004) based and tested this relationship mainly on AGN 
with Log [O m] luminosities of 6.5 and above. Our sam- 
ple has much lower Log [O ill] luminosities, ranging from 
3.8 to 8.3, with median value 5.1. We find that a factor of 
90, instead of 3500, gives a more reasonable value for the 
bolometric luminosity within our sample. This lower fac- 
tor could reflect a higher [O ill] emission relative to L5000 
from our lower luminosity, gas rich systems as compared to 
the SDSS Seyferts. In the case of nuclear radio luminosities, 
we used the sub-arcsecond AGN-related r adio flux to calcu- 
late the im plied jet mechanical power fsee lNagar et al.ll2005l 
for details). iNagar et all l|2005h have demonstrated that this 
method to calculate the bolometric luminosity works well in 
low luminosity AGN. For each nucleus, we then adopted a 
value of Ijtjj^ = Lboi/LEdd (column 8 of Table 3]) from the 
previous four values, keeping in mind their relative uncer- 
tainties and the paradigm that low luminosity AGNs in the 
low (high) state are more likely to produce radio (X-ray) 
emission. This adopted value of the Eddington ratio was 
used to calculate the estimated mass accretion rate (col- 
umn 9). The 'inactive' galaxies NGC 3351, NGC 5248, and 
NGC 5806, all three with HH-type nuclei, have not been de- 
tected in the hard X-ray nor radio. Their Eddington ratios 
are likely lower than the value adopted in Table 3J since a 
significant part of their nuclear [O ill] luminosity could come 
from HII regions. Note that while the overall 'Seyfert' and 



'inactive' samples show expectedly different Eddington ra- 
tios, there is some overlap; two of the Seyferts, NGC 3627 
and NGC 6951, have Eddington ratios similar to the inactive 
nuclei. 

The relationship between mass accretion rate and the 
kinematic of the gas and stars is illustrated in Fig. 1101 The 
result is rather striking: black holes with low accretion rates 
reside in galaxies in which the stellar and gas kinematic axes 
are closely aligned and in which the gas kinematic axis is 
least disturbed, while Black holes with accretion rates higher 
than 10 -4 ' 5 Mq /yr reside in galaxies with either either a 
large misalignment between gas and stellar axes, or a large 
twist in their gas kinematic axes. 

5.2 Feeding the AGN 

Why a fraction of nearby galaxies host an active nu- 
cleus still remains an unsolved question. Several mecha- 
nisms may be involved in the fueling of the central re- 
gions of galaxies, and subsequently activating and sustain- 
ing nuclear activity. Removal of angular momentum and 
transportation of the gas towards the very central region 
may be driven by gravitational perturbations due to in- 
terac tions such as mergers o f gala xies or tidal interac- 
tion (|Barnes fc Hernquistill99ll.ll996h or n on-axisymmetric 
galactic structures i Shlosman et al, |2000| . e.g. bars, spi- 
rals,), as well as turbulence instabilities and shocks in 
the interstellar medium. The majority of optical and near 
infra-red (NIR) imaging studies failed to find any mor- 
phological difference between active and inactive galax- 
ies on a range of scales that encompasse s galactic in- 
terac t ions, stellar bars and nuclear spira ls ( Martini et al.l 
20031; iKnapen et al.l [2000l ; iMarquez et all |2000T ). However. 



Hunt fc Malkan (|2004h found significant excess of isophotal 



twists and morphological disturbances in Seyfert 2 galaxies 
compared to inactive and Seyfert 1 or LINER galaxies. An- 
other recent study of the presence and structures of dust 
in the circumnuclear regions revealed differences between 
early type Seyfert and inactive galaxies l|Simoes Lopes et alj 
2007). These results suggest the possibility of identifiable dy- 
namical differences, in the circumnuclear regions in active 
and inactive galaxies. 

In this study, we investigated the stellar and ionised gas 
kinematics of a well-defined sample of nearby active and in- 
active galaxies, in order to probe the galactic potential in the 
central kpc regions. First, we have measured the kinematic 
misalignments between the stellar and the ionised compo- 
nents. The distribution of kinematic misalignments has of- 
ten been used to co nstrain the or i gin o f the ionised gas i n 
early- type galaxies dBertola et al.lfl992l ; ISarzi et al.l T2006) . 
If the gas is produced internally (e.g. stellar mass loss) and 
comes from the galactic disc, the bulk of the gaseous compo- 
nent in the central region is expected to co-rotate with the 
stars. The angular momenta of the stellar and gaseous com- 
ponents would therefore be globally aligned. Some mild mis- 
alignments may exist due to the effect of non-axisymmetric 
structures (bars, spirals), or the presence of decoupled stellar 
components. Gas acquired from an external source, e.g. via 
interaction with galaxies or the intergalactic medium, should 
eventually settle in the existing principal planes of the grav- 
itational potential. In a triaxial potential hosting a large- 
scale disc, the gas would thus preferentially fall onto either 
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Table 4. Eddington Ratios and Mass Accretion Rates. C olumns are: (1) Pair; (2) Galaxy name; (3) Black Hole mass fro mlLodato fc Bertinl 
i2003T l (NGC 1068, maser dynamics). iKaspi et al] <200C|) (NGC 3227, NGC 4051, NGC 4151; reverberation mapping), fSarzi et al] [|200ll ) 
(NGC 44 59). For the other gala xies, we estimated the black hole mass from the average stellar a in the central 10" of our SAUR0N maps 
using the iTrem aine et all J2002f) relationship; (4) ratio of nuclear (HST) optical luminosity (with no bolometric correction applied; from 
IChiaberge et alj|2005h" to Eddington lu minosity; (5) Eddington ratio derived from the 2-10 keV X-ray luminosity, us i ng th e emperical 
relation found bv lUlvestad fc Hoi <200ll): L b „i = 6.7 Lx (2-10 keV). The data are from the compilation in lPanessa et al" I (120061 s ). except for 
NGC 4459 and NGC 5055 {Gonzalez-Martin et al.ll2006r i and NGC 6951 l|Kim et al.lll992l: from Einstein data with no correction applied 
to transform the luminosity to the 2-10 keV regime); (6) Eddington ratio derived from the [O Hi] luminosity using L bo i = 90 * LjQ m j 
(see text); (7) Ed dington ratio implied by the mechanical energy in the radio jet (Qj e t) as derived from the nuclear radio emission (see 
Nagar et al. 2005]); (8) adopted Eddington ratio: this was roughly calculated from columns 3 to 6 keeping in mind the reliability of the 
various values; (9) the mass accretion rate implied by the adopted Eddington ratio under the standard assumption of a 10% radiative 
efficiency for the accretion process. 



the disc equatorial plane or in a plane perpendicular to it, 
giving rise to a polar ring. Polar ri ngs are observed in a few 
nearby early-type galaxies (see e.g. lBertola et alJll992h . The 
distribution of kinematic misalignments, measured as posi- 
tion angle differences between kinematic major-axes, around 
the reference values of 0, 90 and 180°, should therefore be 
closely linked with the origin of the accreted gas. 

So far, there seems to be no trend for the kine- 
matic misalignments of early-type galaxies with respect 
to Hubble type , gala ctic environment and galaxy luminos- 
ity l|Sarzi et al] 120061 ). In our sample, five Seyfert galaxies 
present significant kinematic misalignments (>20°) however, 
only NGC 2655 presents a gaseous component rotating al- 
most perpendicularly to the stellar component. For all the 
other Seyfert (four objects) and for all the inactive galaxies, 
the ionised gas co-rotates with the stars (kinematic mis- 
alignments smalle r than 20°). This follows the results of 
ISarzi et al] l|2006h for a sample of later-type galaxies and 
therefore suggests that the presently observed ionised gas 
in the circumnuclear regions does not have a purely exter- 
nal origin for both the Seyfert and the inactive galaxies. 
The only clear ex ception is NGC 2655, a galaxy known to 
host a polar ring IjKeel fc Hummellll988l ). It is interacting 
with a companion and had presumably suffered a minor 
merger in the past (|Huchtmeier fc Rich tcr 1982]). Interac- 



tion between galaxies is an efficient way to transport gas 
towards the central parts and to trigger nuclear star for- 
mation l|Barnes fc Hernquistlll99ir i. but it seems neither a 
nec essary nor a sufficient mechanism for most Seyfert activ - 
ity (|Pe Robertis et al.lll998l ; ISchmitdlioOll ; iHo et alj|2003h . 

A more refined analysis revealed that the radial vari- 
ations of the kinematic PAs of the ionised gas are three 
times greater for the Seyfert than for their associated in- 
active galaxies, in the central 1.5 kpc. Indeed, due to the 
relatively large pixel size of SAURON compared with other 
IFUs and the conservative approach of our analysis, the 
intrinsic gaseous deviations may be even larger than sug- 
gested by our results, particularly for the very inner re- 
gions. Although our sample is small, this trend hints for 
the presence of non-axisymmetric perturbations of the grav- 
itational potential in the inner kpc of the Seyfert galax- 
ies. Since the observed stellar velocity fields of these ac- 
tive galaxies are very regular with major axes rather well 
aligned with the outer photometric ones, theses dynamical 
perturbations must be small. The gas, being more sensi- 
tive to small deviations from axisymmetry, responds better 
to small non-axisymmetric structures such as weak bars or 
spirals. This would be in fact in agreement with previous 
imaging studies which found no strong excess of bars or nu- 
clear spirals in Seyfert galaxies. If strong deviations from 
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Figure 10. Distribution of the differences between the global kinematic PAs of the stars and gas (left panel) and the amplitude of 
variations of the gaseous kinematic PAs in the inner 1.5kpc (right panel) as a function of the mass accretion rate (see Table Active 
galaxies are shown as filled red diamonds and inactive galaxies as filled blue circles. For reference, a line is drawn at APA=20° and mass 
accretion rate = 10 — 4 5 Mq/jt in each of the panels. 



axial symmetric were present, the gaseous kinematics would 
have been extremely disturbed and easily detectable with 
standard (long-slit) spectroscopy techniques. The presence 
of weak deviations from axial symmetry in Seyfert galax- 
ies seems to be in good agreeme nt with the refined imaging 
study of iHunt fc Malkanl (|2004j ). They found an excess of 
NIR isophotal twists in Seyfert 2 hosts compared to inac- 
tive galaxies, and discussed the non-axisymmetric instabili- 
ties (nuclear disc, nested misaligned bars, triaxial structures) 
which could be responsible for such morphological pertur- 
bations. This excess of isophotal twist i s not observed for 
the Seyfert 1 galaxies, IHunt fc Malkanl (|2004h interpreted 
this difference between Seyfert 1 and Seyfert 2 galaxies as 
an evolutionary effect, the Seyfert 1 hosts being 'older' than 
the Seyfert 2 ones in their scenario. We would need a larger 
sample to investigate statistical differences in the ionised 
gas properties of Seyfert 1 and 2 galaxies. A larger study 
of Active vs. Inactive galaxies in the SDSS with IMACS is 
now in progress within our group, the active and inactive 
galaxies drawn from the same parent sample so hopefully 
reducing any underlying selection biases and also offering 
more controls per Seyfert (Westoby et al. in preparation). 

In addition to a ready supply of fuel and a transporta- 
tion mechanism, timescales are important. Galaxies are 
thought to experience recurrent episodes of nuclear activity 
in their lifetime, with instabilities in the fueling rate poss i- 
bly regulating the activity cycles l|Saripalli fc Mackll2007h . 
Dynamical instabilities in the inner kpc regions and nuclear 
activity can evolve in parallel. However, due to the time re- 
quired to set up sufficient gas inflow, the non-axisymmetric 
perturbations may be observed while not enough gas has 
been driven in the vicinity of the SMBH and the galaxy is 



inactive. Similarly such gravitational instabilities may have 
disappeared when the nuclear activity is still underway. Ob- 
servationally, it is therefore reasonable to hypothesise that a 
snapshot of active and inactive galaxies might show equal de- 
grees of host galaxy disturbance. However, while disturbed 
gas observed in a host galaxy is not the gas that is simultane- 
ously fueling the nuclear activity, the dynamical time scale 
in the circumnuclear regions is similar to the AGN lifetime. 
Therefore, the correlation between the presence of pertur- 
bation in the ionised gas in the central kpc regions and the 
nuclear activity suggested by our study is striking and sup- 
ports a close link between the AGN, host dynamics on scales 
of hundreds to tens of parcsecs and relevant dynamical and 
activity timescales. 

Finally, the presence of cr-drops in the stellar velocity 
fields observed for some of our gal axies may indicate th e 
presence of cold nucle ar stellar discs (|Emsellem et al.ll200ll ). 
IWozniak etail (|2003T l interpreted such structures as the re- 
sult of star formation following the accretion of gas in the 
central few hundred parsecs. High resolution follow-up IFU 
mapping of the central few arcseconds, using e.g. the OASIS 
IFU, could reveal direct evidence for these discs. Such data 
on a relatively small FOV are rather difficult to interpret. 
However, associated with our larger scale SAURON maps, 
they could extend our understanding of fuel delivery closer 
again to the AGN. 



5.3 Summary 

In this paper, we have presented two dimensional morphol- 
ogy and kinematic maps of the stars and ionised gas compo- 
nent in the central kpc regions of a distance limited sample 
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of well-matched Seyfert and non-Seyfert galaxies. We found 
that: 

• Kinematics of the stellar and gaseous components in 
the circumnuclear regions of our galaxies are dominated by 
disc-like rotation. The stellar kinematics present very regu- 
lar rotation patterns for both inactiv e and Seyfert galaxie s 
in agreement with previous studies IjBarbosa et all 2006). 
Ionised gas is generally co-rotating with respect to the stars. 

• The stellar global kinematic major-axes are aligned 
within 10° with the outer galactic disc major-axis for all 
galaxies, except for two Seyferts: NGC 1068 and NGC 5194, 
for which the determination of the photometric major axis 
orientation is quite uncertain. Except for these two particu- 
lar galaxies, the circumnuclear stellar kinematic major axis 
is a robust determination of the line-of-nodes for either the 
Seyfert and the non-active galaxies. 

• The gaseous component co-rotates with the stars, the 
global stellar and gaseous kinematic PAs being aligned for 
all galaxies within 20°, except for two Seyfert galaxies: 
NGC 2655 and NGC 4579. Our kinematic maps reveal some 
disturbances in the gaseous velocity field, which may be evi- 
dence for deviations from axisymmetry. These perturbations 
in the ionised velocity fields may hint for the presence of 
streaming or radial flows, which could be related to fueling 
mechanisms. 

• The gaseous kinematics PAs have radial variations over 
the SAURON FOV 3 times greater than the stellar ones, for 
either the Seyfert and their associated inactive galaxies. 

• The variations of the stellar and gaseous PAs in the 
central kpc of the Seyfert galaxies are 2.5 times greater than 
those of the inactive galaxies. 

We conclude that the increased incidence of disturbed 
ionised gas in the circumnuclear regions of the Seyferts on 
scales of hundreds to tens of parsecs, compared with inactive 
galaxies, supports a close link between nuclear host dynam- 
ics and relevant dynamical and activity timescales. 
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APPENDIX A: NOTES ON INDIVIDUAL 
GALAXIES 

Al NGC 1068 (M 77) 

NGC 1068 is one of the most nearby Seyfert 2 galaxies 
(D — 14.4 Mpc), and it has been well studied at dif- 



ferent wavelengths. iHo et all (|l997f ) classified this active 
galaxy as a Seyfe rt of Type 1.8, but a dusty torus hides a 
Seyfe rt 1 nucleus l|Antonucci fc MilTerlll985l ; iJourdain et alj 
Il994l ). NGC 1068's n ucleus hosts a radio j et and an ionised 
cone at PA ~ 10° (jMuxlow et al.l 119961 ). ISchinnerer et all 
(2000) summarized the different observed morphological 
structures : the outer disk and the outer oval, the latter 
being interpreted as a primary bar, a HI ring at its Outer 
Lindblad Resonance and the two-arms inner spiral at its In- 
ner Lindblad Resonance. A secondary Near-Infrared (NIR) 
bar extends up to about 16". The SAUR DW data we presente d 
here have been already published bv lEmsellem et ail (|2006h . 
and we refer the reader to this paper for a detailed analy- 
sis: we here only provide a brief description of the revealed 
structures. 

The stellar continuum map (Fig. I4a|l shows elliptical 
isophotes, elongated along the PA of the NIR bar. The stel- 
lar velocity field (Fig. I4ap shows strong departures from 
axisymmetry, with an S-shaped zero velocity line, and a 
slightly varying orientation of the kinematics major-axis. 
The velocity dispersion rises towards the centre, reaching 
200 kms" 1 at R~10"and then presents a d rop in the in- 
ner 5" with values down to about 100 kms -1 . lGerssen et al.l 
(2006) observed this galaxy with the Gemini Multi-Object 
Spectrograph (GMOS) IFU, covering the central 10 x 8 arc- 
sec. Our stellar kinematics properties are in good agreement 
with theirs although they found a kinematic PA offset by 
about 13° from ours, their FOV being too small to detect 
the change of orientation of the kinematic major-axis. 

H/3 and [O ill] distributions are quite different (Fig, l4ap . 
H/3 emission is very high in the inner 5", and traces the spiral 
arms outside this region. As for H/3 , [O ill] emission peaks 
in the central parts, with the distribution of [Oiii] becom- 
ing very asymmetric outwards. It is found predominantly in 
the North-East side, tracing the northern ionisation cone. 
Despite the significant differences between the distribution 
of the H/3 and [O m]line emission, their velocity fields and 
velocity dispersion maps are very similar (Figs. I4"al and lBla} . 
The velocity fields of H/3 and [O ill] both display a prominent 
S-shaped zero velocity curve, evidence for strong deviations 
from circular motions. The velocity dispersion maps of H/3 
and [O ill] show a peak in the central 5" parts and then reach 
lower values outside. As expected, the [Oiii]/H/3 is high in 
the region of ionisation cone (~ 10), and lower in the spiral 
arms ([Oiii]/H/3 < 1), the latter being dominated by star 
formation regions (Figs. l4a)l . 

A2 NGC 2655 

NGC 2655 is an early-type spiral galaxy, hosting a Seyfert 2 
nucleus. Optical and radio spectr oscopic data exhibit of f- 
plane gas in the central part (lErwin fc Sparkel 120031). 
and a dusty polar ring seen in Ha l|Keel fc Hummel 
1988). This galaxy is interacting with its companion 
NGC 2715 and presents remn ant signatures of a past merger 
j|Huchtmeier fc Richterlll982l ). 

The stellar continuum map (Fig. l4b|l shows an elongated 
feature, along the East-West direction, which corresponds 
well with the PA of the line of maximum stellar velocity. 
The stellar velocity field (Fig. I4b[l is remarkably regular, 
with a PA around 90°. The velocity dispersion map shows 
a drop in the inner 3". 
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H/3 and [O ill] emission lines distribution and kinemat- 
ics are similar (Figs. l4bl and lBla)) . The emission lines are very 
bright in the inner 5". Away from the centre, the distribution 
maps show a bright knot to the South-East side, 15" away 
from the nucleus, and a lane at 10" on the West side : both 
features have been reported bv lKeel fc Hummell (|l988l ). The 
SE knot corresponds to high [O iii]/H/3 ratio values (~ 4) 
similar to the ones obsered in th e nuclear region. This fe ature 
could be driven by the radio jet l|Keel fc Hummel 1988). The 
emission lines ratio is smaller in the West lane (around 2- 
3). The emission lines velocity fields show strong departures 
from axisymmetry, with an S-shaped zero velocity curve, the 
kinematic major-axis having a changing PA, from ~ 90° in 
the inner 5" to near 180° in t he outer parts. As hinted by 
Sil 'chenko fc Afanasievl (|2004h . the ionised gas rotates to- 
gether with the stars in the inner 5", but follows the dust 
polar ring at radii > 5" — 6", rotating perpendicularly to the 
galactic plane. [O ill] and H/3 velocity dispersion maps show 
the same features: a rise in the inner 5" (~ 220 kms -1 ) and 
in the western feature (305 kms -1 ), and a roughly constant 
value outside (between 100 and 150 kms" 1 ). 



A3 NGC 2985 

NGC 2985 is an early-type spiral galaxy, with regular gas 
and stellar distributions. It is the control for NGC 4151. 
The disk orientation is determined by an inclination of 
i = 40° and a major-ax is position angle of PA = 180°. 
iNoordermeer et al.l (|2005h report a tidal interaction with its 
companion NGC 3027. 

The stellar velocity field of NGC 2985 shows a regu- 
lar and symmetric pattern, consistent with motions in an 
axisymmetric gravitational potential (Fig. He)) . The PA of 
the kinematic major-axis (-2°) is aligned with the photo- 
metric major-axis of the outer disc. The velocity dispersion 
rises inwards, with 100 kms -1 at 15"from the centre and 
150 kms -1 in the inner 5". There is very little emission 
from the ionised gas H/3 and [O ill] in NGC 2985 (Fig. EeJ|. 
The emission lines flux peaks in the inner 2" with very little 
flux outside. The H/3 (Fig. iBTb] and [O in] (Fig. gej kine- 
matics are similar: the velocity fields are regular, showing 
the same overall symmetry as the stellar velocity field. The 
gaseous velocity dispersion maps show no specific features, 
with values rising towards the centre (200 kms -1 for [Om], 
160 kms -1 for H/3). Finally, the [Om]/H/3 line ratio reaches 
values up to 1.5-2 in the centre and decreases poutwards (for 
r > 10 "). 



A4 NGC 3227 

NGC 3227 is a well studied barred galaxy, interacting with 
the elliptical galaxy NGC 3326. and with an associated 
dwarf galaxy l|Mundell et al.| [ 20 04) . The galactic disc has an 
inclination of 56° , with a o uter photometric ma jor-axis at a 
position angle PA = 158° l|Mundell et al.lll995h . coincident 
with the major- axis of the stell ar bar. It is hosting a type 1.5 
Seyfert nucleus (|Ho et al.lll997h . The centr al region has been 
mapp ed in 12 CO(1-0) and 12 CO(2-l) by ISchinnerer et all 
(2000), who detect molecular gas very close to the nucleus 
(~ 13 pc). The inner kpc region s host several comple x fea- 
tures: a radio jet (PA ~ -10°; iMundell et~ai1 (|l995r )1, an 



conic al NLR outflow at a PA of about 15° (Mundc ll et al 



1995), an Ha outflow at PA ~ 50 ° dArribas fc Mediavilla 



1994) , and a molecular nuclear ring ( Sch innerer et al.ll2000l V 
Our SAURON stellar continuum map (Fig. I4a[) presents 
an elongated structure (PA = 153°), aligned with the stellar 
bar. The stellar velocity field is regular (Fig. I4ap . the major- 
axis being parallel to the galactic disc orientation. Outside 
r > 5", the stars rotate rigidly (isovelocity contours parallel 
to each others). The stellar velocity dispersion rises towards 
the centre, reaching 200 kms -1 in the inner 5". 

Fig. [4a] presents the [O ill] and the narrow component 
of H/3 emission lines. The central broad component for the 
H/3 line corresponding to nuclear emission (BLR) has been 
removed. [O ill] distribution is rather asymmetric, elongated 
on the North side of the field and H/3 emission lines extends 
from South-East to East. The [O ill] and H/3 velocity fields 
are similar: they show strong deviations from axisymmetry 
with a PA of the major-axis changing from ~ 30° in the 
central regions, to about 170°, close to the stellar kinematic 
major-axis orientation. The velocity dispersion of both H/3 
and [O ill] rise towards the nucleus, reaching 310 kms" 1 and 
360 kms -1 , respectively in the centre. Finally, the [O m]/H/3 
lines ratio map shows an elongated structure of high values 
(> 7) along the global kinematic major-axis (North-South). 

A5 NGC 3351 (M 95) 

NGC 3351 is a barred inactive spiral galaxy. It is the con- 
trol object for NGC 4579. NGC 3351 hosts a large-scale stel- 
lar bar and an inner mol ecular bar-like feature in the cen- 
tre l|Devereux et al.lll992f ). The latter molecular structure is 
aligned with the major-axis of the galaxy (PA ~ 10°) and 
perpendicular to the large scale bar. Two rings of HII regions 
associated with resonances due to the stellar bar exist: one 
inner ring with a radius of 10" and anoth er at 70" from the 
centre. More recent observations in CO (|Joeee et al.l l2005h 
show evidence of the non-circular streaming motions in the 
inner region. 

The SAURON stellar continuum map (Fig. [4f|| shows a 
ring- like structure with 6" < R < 10", which corresponds 
well with the HII ring at the presumed ILR of the large scale 
bar. The stellar velocity field is regular throughout the field 
of view. Stellar velocity dispersions are slightly lower (by 
~ 20 kms -1 ) inside the ring than outside. 

[O ill] and H/3 distributions and kinematics do not differ 
much from each others. The intensity maps clearly show 
the known HII ring at ~ 10" (Fig. [4f]), and there is very 
little emission outside that structure. The gas velocity maps 
are regular, following the orientation observed in the stellar 
velocity field. The [O ill] velocity dispersion map is almost 
featureless, with only a slight gradient from North-East (100 
kms -1 ) to South-West (140 kms -1 ). There is a velocity 
dispersion drop (~ 50 kms -1 ) for H/3 in the ring (Fig. lBlb)) . 
The [O ill] /H/3 lines ratio map show low values in the inner 
regions (< 0.1), as expected from star formation. [Om]/H/3 
ratio is significantly higher outside the ring (~ 1), probably 
mostly due to the uncertainty in the emission line flux there. 

A6 NGC 3627 (M66) 

NGC 3627 is an SAB spiral galaxy hosting a Seyfert 2 nu- 
cleus. It is part of the Leo Triplet system, interacting with 
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NGC 3623 and with evidence for a past interaction with 
NGC3628 (|Zhang et alj|l993h . Previous radio and optical 
observations show asymmetric structures in the velocity 
field, non -circular gaseous mo tions and a molecular CO in- 
ner ring (IChemin et alj|2003h. Mo r e rece ntly, IFS data ob- 
tained by lAfanasiev fc Sil'chenkd (|2005f ) showed evidence 



for strong gaseous radial motions. 

The SAURON stellar continuum map shows an elongated 
structure along the photometric PA ~ 170° (Fig. He)) . The 
velocity field presents a regular rotation pattern, with a 
slight twist of the zero-velocity line in the very central parts 
(r < 3"). The stellar rotation follows the orientation of the 
photometric major-axis. The stellar velocity dispersion rises 
towards the central regions up to 120 kms -1 . 

The [O ill] emission is concentrated in the inner 5", with 
some emission in an elongated structure, along the photo- 
metric PA (Fig. I4c[) . As for [Om], H/3 emission is prepon- 
derantly observed in the central parts, with an additional 
bright spot about 15" South of the nucleus. The [O ill] 
and H/3 kinematics are quite similar (Figs. [4c] and IBla|) . 
Their velocity fields show strong deviations from axisym- 
metry: the orientation of the kinematic major-axis changes 
from the external regions where it is aligned with the stellar 
kinematics, toward the nucleus where the ionised gas kine- 
matic major-axis devia t es by ~ 40° from the stellar one. 
lAfanasiev fc Sil'chenkol (|2005l ) suggested that strong radial 
motions may exist in the central parts of this galaxy. Both 
[O ill] and H/3 velocity dispersion maps are almost feature- 
less, with a slight depression in the very centre. Finally, the 
[O ill]] /H/3 lines ratio reaches its highest values in the inner 
2". 



A7 NGC 4051 



NGC 4051 is an SAB galaxy classified as SI. 2 by iHo et all 
{ 19971 ). The stellar bar extends up to ~ 50" along a PA 
ss 1 35°. Radio-observations sh ow a triple source at a PA ~ 
73° (|Christopoulou et al.lll997T ). and the [Om] emission line 
distribution in the inner 3" is aligned with this radio com- 
ponent. The narrow-line profiles of NGC 4051 present strong 
blue wings. This blue-shifted component is also detected 
with SAURON in the H/3 a nd [Q hi] lines, and has th erefore 
been fitted (see Sect. l2.4p . lChristopoulou et all (|l997l "t found 
evidence for a conical outflow, at l'.'5 from the nucleus. 

Our SAURON stellar continuum map (|4d[) show regular 
isophotes elongated along the PA of the large-scale bar. The 
stellar velocity field presents a regular and symmetric rota- 
tion pattern, the rotation axis being aligned with the photo- 
metric major-axis (PA~ 135°). The stellar velocity disper- 
sion decreases in the inner 5"down to ~ 50 kms - . 

The emission line intensity is slightly extended towards 
the North-East, c onsistent with the io n ised [ O hi] outflow 
cone observed by IChristopoulou et al.l (|l997 ). The orien- 
tation of the ionised gas kinematic major-axis varies with 
radius: it is aligned with the stellar one in the outer part 
of the SAURON FOV, and then deviates by about 30° from 
the line-of-nodes of the galaxy. The H/3 emission line ve- 
locity dispersion reaches its highest values in the inner 5", 
while a decrease of dispersion is observed in these central 
regions for [Om]. The [O ill] /H/3 lines ratio peaks at about 
3" North-East from the nu cleus, within the outflow region 
l|Christopoulou et al.lll997h . 



A8 NGC 4151 

NGC 4151 is one of the most well-studied galaxies of our 
sample. For a review on this galaxy and the AGN proper- 
ties, see lUlrichl <|2000t ). It is an almost face-on bar red spiral 
galax y hosting a Seyfert type 1.5 active nucleus l|Ho et al.l 
Il997h . The large-scale weak bar is elongated along a PA 
of about 130° and the photometric m ajor-axis orientation 
is close to 20° (|Mundell fc Shone! fl999?). Its nuclear contin- 
uum and BLR emission are highly variable. It contains a 
radio jet along a PA of 77° and its Narrow Line Region 

NLR) is extending over ~ 10", and is aligned with the jet 

Mundell et al.ll2003l ). 

Our stellar continuum map (Fig. lie)) presents rather 
round isophotes. The stellar velocity field shows a twist of 
the zero-velocity line in the inner 5 " which corresponds to 
the region where the BLR dominates the spectral features. 
A er-drop is observed in the central part of the velocity dis- 
persion map, but could be due to the BLR contamination. 
The [O ill] and H/3 emission lines distributions are elongated 
on the South- West direction at a PA of about 50° associ- 
ated with the N LR and the Extended NLR (ENLR, R > 5" 
lAsif et al.|[2005l ). The emission lines ratio reaches very high 
values (~ 12) in this region. The ionised gas kinematics is 
dominated by rotation, but with some disturbances on the 
North- West side. [O ill] emission line velocity dispersions 
peak in the inner few arcsec (~ 280 kms -1 ). Outside the 
very central parts, the velocity dispersion is between 50 and 
130 kms -1 . It increases at the edge s of our map, correspond - 
ing to the locati on of a dusty ring l|Vila-Vilaro et allll995l ). 
lAsif et al] l)l998[ ) revealed the presence of streaming motions 
along this dust structure, which may be associated with the 
high [O ill] velocity dispersion observed in our SAURON maps. 
H/3 velocity dispersion is flat in the outer parts (R > 10", 
with a between 80 and 130 kms -1 ) and rises regularly to- 
wards the center, reaching 200 kms -1 in the inner 3". The 
maximum velocity dispersion (~ 350 kms -1 ) is reached at 
about 3" 5 North-East of the centre. 



A9 NGC 4459 



This SO galaxy is in the Virgo cluster. It is the con- 
trol non-active galaxy for NGC 2655. The SAURON data 
resented here (Fig. I4b|l are part of the SAURON project 
de Zeeuw et al j|2002h and were therefore already published 
elsewhere (|Emsellem et~ai1l2004l ; ISarzi et alj|2006j ). 

The stellar component of NGC 4459 shows roundish 
isophotes, and its stellar velocity field presents a very reg- 
ular rotation pattern the major-axis of which is parallel to 
the line-of-nodes (PA ~ 77°). The velocity dispersion rises 
towards the centre, reaching 200 kms -1 (Fig. I4b() . There is 
almost no ionised gas emission (our maps have been clipped 
at 3a), except for the very centre (R < 3"). The emission 
lines velocity field seems to be dominated by rotation, and 
the velocity dispersion rises from about 80 kms -1 for R > 3" 
to 160 kms -1 for H/3 (200 kms -1 for [Om]) inside 3". The 
[O iii]/H/3 line ratio rises regularly towards the centre. 



A10 NGC 4579 (M 58) 

NGC 4579 is an active galaxy classified as S1.9/L1.9 
l|Ho et al]|l997h . It hosts a NIR stellar bar of 9 kpc of diam- 
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eter, oriented along PA = 58° . iGarcfa-Burillo et alj (|2005l ) 
observed this galaxy with the Plateau de Bure Interferome- 
ter (CO) as part of the NUGA survey. They found a nuclear 
molecular spiral from R ~ 1 kpc down to ~ 200 pc, driven 
by the stellar bar and detected highly non-circular motions 
over the spiral arms. They interpreted these perturbations 
as outflow motions. 

Fig. |4f| presents our SAURON maps for this galaxy. The 
stellar continuum isophotes are rather round in the centre 
and slightly elongated along a PA of ~ 60° corresponding to 
the orientation of the stellar bar. The stellar velocity field 
presents a very regular rotation pattern, with a kinematic 
major-axis oriented at a PA of 95°. The stellar velocity dis- 
persion rises regularly towards the centre. 

[O ill] and H/3 distribution and kinematics are very 
similar (Figs. [4f] and IBlb|l . The ionised gas intensity maps 
show a spiral-like structure which cor responds to the nu- 
clear molecular structure observed by iGarcia-Burillo et al.l 
(2005). The [Om]/H/3 ratio presents lower values in the spi- 
ral arms than outside. The ionised gas kinematics is com- 
plex and shows strong departures from axisymmetry. The 
gas kinematic major-axis (PA — 100°) is almost perpendic- 
ular to the stellar bar, and varies within the SAURON field of 
view. Despite the similarities between H/3 and [O ill] distri- 
butions and velocity fields, their velocity dispersions are sig- 
nificantly different. The U[oiii] map presents an elongated 
structure of high values along the stellar bar (Fig. [4f| , while 
°Fl/3 decreases in the central parts Fig. IBlbl 



All NGC5055 (M63) 

NGC 5055 is a spiral non-active galaxy, control for 
NGC 5194. Previous radio and optical studies show an over- 
all remarkable r egularity and symmetry in its morphology 
and kinematics ([Thornlev fc Mundvl 1 19971 : iBattaglia et al.l 
2006). There are evidence for a tidal interaction with a 
compani on (UGC8313), and the out er disk of NGC 5055 is 
warped. iBlais-Ouellette et al.l l|2004h studied the Ha emis- 
sion of NGC 5055 with a Fabry- Perot: they found two veloc- 
ity components in the inner 8"of this galaxy. One component 
is consistent with the global disk kinematics while the sec- 
ond one exhibits a pattern compatible with a bipolar outflow 
or a counter-rotating disk. 

Our SAURON stellar continuum map presents an asym- 
metry towards the North-West (Fig. |4g[ ). The stellar velocity 
field shows very regular rotation pattern with the kinematic 
major-axis aligned with the photometric one at a PA of 100° . 
The velocity dispersion rises in the very centre (r < 2"). 

The ionised gas H/3 and [O ill] are mainly concentrated 
in the inner 5". The [Om] distribution shows a V-shaped 
structure towards the North, while the H/3 intensity map 
shows a number of hot spots broadly distributed over the 
FOV. The [O ill] and H/3 velocity fields are very regular 
and similar to the stellar velocity field (except for their 
higher amplitudes) . [O in] velocity dispersion maps is fea- 
tureless but H/3 velocity dispersion map shows an elongated 
structure of high dispersion roughly aligned with the kine- 
matic major-axis. Finally, the [Om]/H/3 lines ratio reaches 
its highest values in the region of high [O ill] emission line 
flux. 



A12 NGC 5194 (M 51) 

NGC 5194 is an almost face-on spiral galaxy in interaction 
with NGC 5195 . Its nucleus is classified as a Seyfert type 2 
l|Ho et al.| [l9971. Many studies have been made on the struc- 
ture and the nucleus of this nearby galaxy. A large outflow 
bubble extents up to ~ 9" (~ 42 pc) North of the nu- 
cleus and a bright cloud is located about 3"South of the 
nucleus. The se structures hav e been observed at different 
wavelengths (|Ford et alj [l985). and they are found to trace 
the biconical NLR along a PA of 163°. NGC 5194 is hosting 
a radio-jet which interacts with the surrounding interstellar 
medium. 

The SAURON stellar continuum map (Fig. |4f| shows the 
presence of dust lanes. The stellar velocity field has a regular 
rotation pattern with its major-axis at a PA ~ —165°, devi- 
ating by ~ 2 5° from the presumed line of nodes orientations 
(PA = 170°. lTullvlll974h . The stellar velocity dispersion is 
featureless. 

H/3 emission line flux is very high in the centre and for 
R > 5" it traces the spiral arms of the galaxy. The [O ill] 
emission line flux traces the ionisation cone of NGC 5194. 
The [O ill] intensity map shows the outflow bubble in the 
North and the bright cloud in the South of the nucleus cor- 
responding to the end of the radio-jet. The emission line 
ration [O m]/H/3 peaks in the centre (with values up to 8) 
and in the southern cloud (values between 3 and 5). The 
ionised gas velocity field seems to be dominated by simple 
galactic rotation, but presents some distortions very proba- 
bly due to the outflow observed in the inner 10". The ionised 
gas velocity dispersion rises towards the centre. 



A13 NGC 5248 

NGC 5248 is a nearby inactive spiral galaxy, control for 
both NGC 4051 and NGC 6951. It har bours two prominent 
circum-nuclear rings l|Laine et al.ll200ll ) and an inner molec- 
ular spiral within 1.5 kpc driven by the large scale bar 
ijjogee et alj|2002h . The a rms of the sp i ral co nnect with the 
inner circumnuclear ring. iJogee et all (|2002h also mention 
the presence of streaming motions in this galaxy. 

The stellar continuum shows an elongated structure 
along the global disk orientation at PA ~ 150° (Fig. I4d[) . 
The velocity field presents a regular rotation pattern, the 
major-axis of which is oriented at a PA of 1 15°, aligned 
with the photometric major-axis (PA = 105 deg lJogee et al.l 
2002). A a-drop is revealed by our velocity dispersion map 
in the inner 3" with an amplitude of ~ 60 kms -1 . 

The H/3 emission line distribution shows a ring-like 
structure with a radius of ~ 5". The [Om] distribution 
is concentrated within the central few arcseconds, and in 
the two brightest peaks of the ring. This ring-like struc- 
tu re corresponds to t he known ring of HII region described 
bv lLaine et al.l (|200ll ). The [Oiii]/H/3 emission line ratio is 
lower inside the ring (0.15), consistently with the presence 
of star formation, and rises slightly towards the very centre 
(0.7). The ionised gas velocity field is dominated by rotation 
with its major-axis aligned with the stellar kinematic one. 
The H/3 velocity dispersion map shows a depression in the 
star formation ring (50 kms -1 ), and slightly rises towards 
the very centre (80 kms -1 , Fig. IBlb|) . The [Om] veloc- 
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APPENDIX B: H/3 KINEMATIC MAPS 

In this section we present the H/3 kinematic maps for each 
galaxy. 

This non- act iv e gal axy is the control for NGC 3627. 
ICarollo et al.l (|2002T ) observed this galaxy with the 
HST/NICMOS2 camera in the near-infrared. They detect 
a prominent bulge and a stellar bar structure in the cir- 
cumnuclear regions. They classified NGC 5806 as a galaxy 
with a concentrated nuclear star formation mixed with dust. 
Our SAURON data is presented in Figs, gc] IBlbl and . The 
stellar morphology and kinematics are quite regular. The 
stellar continuum within the SAURON field of view is elon- 
gated along the outer photometric major-axis of the galaxy 
(PA = 170°). The stellar velocity field shows a regular ro- 
tation pattern with the kinematic-major axis aligned with 
the photometric major-axis. The velocity dispersion is rising 
toward the nucleus reaching about 120 kms -1 in the centre. 

The [O ill] and H/3 emission line distribution show elon- 
gated features along the photometric PA. That feature ex- 
hibits a ring-like structure in H/3, broken in the North- East 
part, at a radius ~ 3". The [O iii]/H/3 ratio shows rather 
low values in the H/3ring (~ 0.3), consistent with star for- 
mation. The gas velocity maps are dominated by rotation, 
and resemble the stellar velocity field, though with a higher 
amplitude. The ionised gas velocity dispersion maps show 
no particular structures. 

A15 NGC 6951 

NGC 6951 is a spiral galaxy which has a large s tellar bar 
(PA ~ 90° ). Its nucleus is cl assified as Seyfert 2 (|Ho et al.l 
1 19971 ). but iPerez et al.l i|2000r ) suggest that it could be con- 
sidered as a transition object, between a very high excita- 
tion LINER and a Seyfert 2. The ionised gas (Ha) is mainly 
concentrated in the nuclear region and form an annulus of 
radius ~ 5" , in wh ich massive star formation is taking place. 
IPerez et alj (|2000T ) interpreted the Ha kinematics in the cir- 
cumnuclear zone as a series of nested disks, decoupled from 
one another. 

The stellar continuum presents isophotes elongated 
along a PA consistent with t he galactic disk major-axis 
iRozas et all (PA- 138° 120021 ). The stellar velocity field 
shows regular rotation pattern whit the major-axis aligned 
with the line-of-nodes of the galaxy. The stellar velocity dis- 
persion shows a drop in the inner 5". 

The H/3 emission line distribution map shows a inner 
broken ring elongated along the major-axis of the galaxy 
(PA = 137°) while [Om] is dominated by the AGN contri- 
bution, its intensity maps showing a compact central peak 
in the inner 3", with no particular feature outside. The emis- 
sion line ratio [O in] /H/3 reaches very low values in the ring 
(0.2), consistent with the presence of on-going star formation 
in the annulus, and rises towards the centre (values of up to 
6) where the effect of the AGN is prominent. The velocity 
fields are quite regular, dominated by the global rotation, 
and aligned with the stellar kinematics. The emission line 
velocity dispersion rises rapidly towards the centre reaching 
about 100 kms -1 close to the nucleus for both the H/3 and 
[O ill] emission lines . These values are in good agreement 
with those found bv lPerez et"al] (|2000l ). 



ity dispersion map shows not specific features, with roughly 
constant values around 50 kms -1 . 

A14 NGC 5806 
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Figure Bla. Velocity and velocity dispersion of H/3 emission lines. Each row present a pair : Seyfert on the left, Control on the right, 
except for NGC 1068 and NGC 3227 which are shown together in the first row. For NGC 4459, the [Oiii] velocity and velocity dispersion 
maps are shown (see Sec. [3] All the maps are orientated so that the photometric major-axis of the galaxy is on the x-axis. The units are 
in kms -1 . The color scale is shown on the right hand side and the cut levels are indicated in the top right corner of each map. The long 
and short arrows on the right of each galaxy name show the North and the East directions, respectively. 
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Figure Bib. Fig. IBlal continued 
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